The origin of stiffening in cross-linked semiflexible networks 
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Strain stiffening of protein networks is explored by means of a finite strain analysis of a two-dimensional 
network model of cross-linked semiflexible filaments. The results show that stiffening is caused by non-affine 
network rearrangements that govern a transition from a bending dominated response at small strains to a stretch- 
ing dominated response at large strains. Thermally-induced filament undulations only have a minor effect; they 
merely postpone the transition. 
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There is a deep interest in the mechanical response of bio- 
logical tissues and gels in view of the importance for biolog- 
ical functions such as cell motility and mechanotransduction. 
Many network-like biological tissues respond to deformation 
by exhibiting an increasing stiffness, i.e. ratio between change 
of stress and change of strain. This has been demonstrated by 
micropipette and microtwisting experiments 1 1 ] on individ- 
ual cells and through rheological experiments on in-vitro gels 
of cytoskeletal filaments (actin, vimentin, keratin |2, 3, 4] and 
neuronal intermediate filaments |5]), as well as on fibrin OLL- 
These biological gels fall within the class of semiflexible poly- 
mers, which has also attracted much theoretical attention in 
the last decade However, these theoreti- 

cal studies have primarily focused on the small-strain regime, 
tractable for analytical treatment. 

In a simple conceptual view, a biopolymer network is an 
interlinked structure of filaments. Thus, stiffening can re- 
sult from stiffening of the polymeric filaments between cross- 
links, from alterations in the network structure, or both. The 
current paradigm is that stiffening is primarily due to the stiff- 
ening of the filaments themselves. This idea has been worked 
out in detail very recently by Storm et al. Q] by adopting 
the worm-like chain model for actin filaments in combina- 
tion with the assumption that the network deforms in an affine 
manner, i.e., each filament is assumed to follow the overall de- 
formation. The worm-like chain model is a well-documented 
description for the stretching of semiflexible polymers, where 
the longitudinal stiffness of undulated filaments is attributed 
primarily to bending; the axial stiffness of the polymeric chain 
itself is much higher am. As the filament is stretched (at 
constant temperature), the amplitude of the transverse thermal 
undulations reduces and, as a consequence, the stiffness in- 
creases. In the limit that the filament is pulled straight, all sub- 
sequent axial deformation would have to originate from axial 
straining of the chain, but at an enormous energy cost. Given 
this description of individual filaments, Storm et al. pro- 
ceed by considering a network consisting of infinitely many 
filaments. Initially the filaments are randomly orientated, and 
as the sample is deformed the network is assumed to distort in 
an affine manner. The affine deformation assumption is well 
known in network models for rubber elasticity, and allows for 
a relatively simple description of the overall network response 
on the basis of the behavior of a single filament. The small- 



strain affine deformation assumption in two-dimensional net- 
works of straight filaments has recently been studied in great 
detail by Head et al. II ill , who conclude that its validity de- 
pends on the cross-link density of the filaments. Their conclu- 
sion, however, cannot be immediately transferred to the stiff- 
ening results of Storm et al. since it applies only to the initial 
response. 

In contrast with the cited literature, we demonstrate in 
this paper that stiffening lies in the network rather than in 
its constituents. During deformation, the filaments rotate in 
the direction of straining, which induces a transition from 
a bending-dominated response to one that is controlled by 
stretching of aligned filaments. By comparing cross-linked 
networks with and without thermal undulations, we show that 
filament reorientation is the dominant mechanism, while the 
thermal undulations only postpone the onset of stiffening. 

Our model is a two-dimensional network model of fila- 
ments in a periodic unit cell of dimensions ffxW. The net- 
work is generated by randomly placing filaments with length 
L at random orientations inside the cell, with proper account 
of periodicity. Thermal undulations are mimicked by super- 
posing on each filament transverse normal modes of the type 
b n sin(nKx/L), where the ampitudes b„ follow a Gaussian dis- 
tribution (cf. l3) with standard deviation y/2/{l v L) (L/nn) 2 . 
The persistence length l p , i.e., the distance over which the fil- 
aments appear straight, is expressed as l p oc k/ (k^T) in terms 
of the bending stiffness K of the filaments, Boltzmann's con- 
stant &b and temperature T. Clearly, the filaments are straight 
in case the temperature is low or the bending stiffness is high. 
We use the first 10 normal modes to generate the initial ge- 
ometry of the filaments and treat / p as an independent quan- 
tity. During mechanical loading, thermal effects are no longer 
taken into account. 

Points where filaments overlap are considered to be cross- 
links, similar to thejprocedure used by Head et al. flUl and 
Wilhelm and Frey II 211 . The networks generated by this pro- 
cedure are taken as the initial, stress-free configuration. In 
the calculations, the cross-links are assumed to be stiff, so 
that both displacement and rotation of the two filaments at 
the cross-link point remain the same. The filaments are elastic 
rods, characterized by a stretching stiffness fj (axial force fl^l 
needed to induce a unit axial strain) and a bending stiffness K 
(bending moment needed to induce a unit radius of curvature). 
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For isotropic elastic rods these values are related through their 
cross-sectional geometry, but are treated here as independent. 
The density of the network is characterized by the line den- 
sity p, i.e., the total length of filaments in the unit-cell divided 
by the cell area, W 2 . For networks with straight filaments 
(Zp/L — * oo), the average distance between cross-links, Z c , is 
inversely proportional to p through Z c — n/p | 16J. We con- 
sider networks of different densities, but only above the rigid- 
ity percolation threshold 1 12]. 

For the numerical study we use the finite element 
method, discretizing each filament with 10 equal-sized Euler- 
Bernoulli beam elements accounting for stretching and bend- 
ing. Geometry changes are accounted for by an updated La- 
grangian finite strain formulation. All filaments are perfectly 
bonded to rigid top and bottom plates, with the top plate dis- 
placed horizontally relative to the bottom plate over a distance 
FW, corresponding to an applied shear strain F. The macro- 
scopic shear stress x is calculated from the total horizontal 
reaction force at the top, divided by W. Convergence studies 
ensured that the cell size W does not affect the results. 

The parameters governing the system are x, T, /j, K, L, p 
and Z p . We choose to present the results through the follow- 
ing dimensionless parameters: x = tL//j, F, p = pL, Zp = l p /L 
and Z~b = \/ K / (l>L 2 ). Note that Zj, is a measure for the ' floppy - 
ness' of the filaments, which reduces to the slenderness ratio 
(thickness over length) for isotropic elastic rods, and / p sets 
the initial shape of the filaments. 




FIG. 1: Average stress (x) versus strain (r) response for networks 
with 4 = 2.3 x 1CT 4 and jiz/L = 1.6 MPa at three different densities 
p and for straight (Z p >> 1) and undulated (Z p = 1) filaments. The 
error bars have a length of two times the standard deviation in strain 
at a given T for ten different realizations at each density. The squares 
correspond to three instances during deformation for which the net- 
work geometry is shown in Fig.[2| 

In a first set of calculations we take Zj, = 2.3 x 10~ 4 and 
fi/L =1.6 MPa (representative for actin microfilaments Jnll ') 



with a density of p = 13, which is well above the rigidity per- 
colation threshold of p = 5.7 [12]. The persistence length is 
taken to be much larger than the filament length (Zp >> 1), 
corresponding to straight filaments. Figure^shows the stress- 
strain response (averaged over ten different random realiza- 
tions). Three regimes can be identified: a regime with a rel- 
atively low stiffness dx/dF, a transition regime and a high- 
stiffness regime. Figure|2]shows three snapshots of the p = 13 
network geometry at F = 0, 0.08 and 0.24 for a typical real- 
ization close to the average response (see the solid squares in 
Fig-[D- Comparison of Fig.|2j) with Fig. |2^ reveals that many 
initially straight filaments have deformed by bending, which 
corresponds to the characteristic low stiffness at small strain 
levels for these densities II ll ll 211 . Subsequently, during the 
transition regime, percolations of stretched filaments appear 
that connect the top and bottom of the cell along a ~ 45° di- 
rection, Fig.|2j;. These filaments are loaded in axial tension, 
resulting in a higher overall stiffness. Thus, Figs. |2j) to |2j; 
demonstrate the transition from a bending dominated regime 
at small strains (total mechanical energy of the systems pri- 
marily consists of bending energy) to a stretching dominated 
regime at higher strain levels (total energy dominated by the 
axial stretching energy). 

The stress-strain response for two higher densities, p = 19 
and 25, is included in Fig. [0 The standard deviation in 
strain from ten realizations is approximately independent of 
the stress value for a given density. The scatter in strain, de- 
fined as the ratio between the standard deviation and the av- 
erage, is independent of the density and has a value of ap- 
proximately 0.2. Figure [0 shows that a certain stress level 
is achieved at smaller strains in case the network is denser. 
The network thus gets stiffer with increasing density, while 
the transition from bending to stretching becomes less abrupt 
and occurs at smaller strain levels. 

Next, the effect of thermally-induced undulations is inves- 
tigated. The same parameters are used as before, but now we 
use Zp = 1, in accordance with experimental findings for actin 
filaments 11711 . Changing the persistence length from Z p ^> L to 
Z p = L is physically similar to increasing the temperature from 
OK to 293K before cross-linking and loading the network in- 
stantaneously. Figure [3] depicts the initial geometry for a net- 
work where the initial filaments' end-to-end vectors have the 
same location and orientation as shown in Fig.|2^. The stress- 
strain results included in Fig. ^ (dashed lines) show that the 
thermal undulations do not change the shape of the overall 
stress-strain curve, but merely delay the transition from bend- 
ing to stretching. The associated 'delay' strain at an applied 
strain level of F = 0.25 is 0.018, 0.029 and 0.024 for p = 13, 
19 and 25, respectively. Similar values were found for densi- 
ties up to 38. Clearly, only a small fraction of the total strain 
is due to the presence of thermally-induced undulations. 

To make connection with the small-strain study by Head 
et al. II 1H . we monitor the degree of affinity of the network 
during straining to large deformations. For this purpose, we 



(a) (b) (c) 

FIG. 2: (a) Initial, r = 0; (b) intermediate, T = 0.08 and (c) large strain, T = 0.24, network configurations for a typical p = 13 realization close 
to the average response shown in Fig.Q(squares). 




FIG. 3: Initial conformation of a p 
filaments corresponding to / p = 1 . 
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define the deviation from affine behavior, AA, as 
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where ||r|| 2 = r ■ r, n the number of cross-links and rW is the 
current position vector of cross-link k. Ar™ is the increment 
in the position of cross-link k during a shear increment Ar in 

(k) 

the simulations, while Ar^ is the corresponding value were 
the deformation affine. Figure@]shows the evolution of AA as 
a function ofT for the cases shown in Fig.^ It is observed that 
the deformation is not affine at small strains, in accordance 
with | ljj, ll 211 . while the deformation becomes increasingly 
affine (AA — > 0) with increasing strain. By comparing Fig.|4] 
with Fig.^ we find that in the transition from the bending to 
the stretching regime, AA increases significantly, indicating a 
reorientation of the filaments. The peak in the AA~r curve 
occurs for both straight and undulated filaments in the transi- 
tion regime. This is another indication of the fact that filament 



undulations do not change the nature of network deformation, 
but merely enhance the strain value at which stretching sets 
in. Once stretching has set in, the deformation becomes more 
and more affine. 
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FIG. 4: The deviation from affine behavior, AA, as a function of 
strain, T, for ij, = 2.3 x 10~ 4 and /u/L= 1.6 MPa, at three different 
densities and for straight (/ p >> 1) and undulated (/ p = 1) filaments. 

To study the influence of the filament properties, the calcu- 
lations are repeated, but with a larger bending and stretching 
stiffness, — 8.3 x 10~ 4 and/j/L = 8 MPa (representative for 
microtubuli lllilo . Note that the bending stiffness increases by 
a factor of 65, while the stretching stiffness becomes five times 
larger. Because of the enhanced bending stiffness, the per- 
sistence length at the high temperature increases from l p = 1 
to 65 so that the filaments are almost straight (in accordance 
with experimental observations on microtubuli 1 17]). Figure[5] 
shows the instantaneous shear stiffness, G = dz/dF, as a func- 
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However, due to non-linear geometrical effects at large strains, 
the stiffness increases with straining according to 1 20] 

G°cpr 2 (i + rr 3/2 . (2) 

This scaling is seen from Fig. [3] to successfully characterize 
the steady state stretching regime at large strains. 

This study leads to the following conclusions. Stiffening 
of cross-linked semifiexible networks is caused by the tran- 
sition of a bending-dominated response at small strains to a 
stretching-dominated response at large strains. This transition 
is mediated by network rearrangements that are not affine. 
Filament undulations only have a minor effect; they merely 
postpone the transition from bending to stretching. Above 
a density-dependent transition strain, the network stiffness 
scales linearly with density and the filament's stretching stiff- 
ness, which is expected to hold in three dimensions as well. 



FIG. 5: Shear stiffness as a function of strain for networks with 
straight filaments having 1\, = 2.3 x 10 or 8.3 x 10~ 4 . The dot- 
ted lines represent the scaling relationship J2j, fitted to p = 25. 

tion of shear strain T for the two different values of % (and ref- 
erence stress fj/L). It can be observed that for the 'floppy' fil- 
aments (7b = 2.3 x 10~ 4 ) the transition from the low-stiffness 
to the high-stiffness regime shifts to lower strains with in- 
creasing density, consistent with Figs.^and[4] For densities 
higher than p = 25, the transition from bending to stretching is 
no longer accompanied by severe filament reorientations that 
cause a peak in AA (see Fig.|4j, but progresses more smoothly. 
At higher densities of either floppy or stiff filaments, the stiff- 
ness first decreases with strain at small strain levels. This is 
caused by buckling of filaments oriented at 135° away from 
the horizontal axis (positive to the right in Fig. [3, which are 
loaded primarily in compression. Figure [5] also shows that at 
small strains the overall stiffness for the floppy filaments is 
much lower than that of the stiffer filaments, but converges to 
the same value at larger strains. This reflects that bending is 
the dominant deformation mode in the low-stiffness regime at 
small strains, while filament stretching governed by /j domi- 
nates at large strains. 

When filament stretching is the dominant deformation 
mechanism, the following scaling arguments hold. The over- 
all stress x is distributed, on average, over the filaments 
through the force F °< xl c , where l c « 1/p is the 'mesh size' of 
the network. This force causes the filaments to deform axially, 
yielding an elongation 8 °c Fl c /fJ, which yields an average 
strain r <* 8/Z c - Substitution yields r«x/ (p/j) so that G <* p/j. 
In three dimensions the same scaling relation holds 1 19[. It 
should be noted that in the high-stiffness regime, the network 
'locks' when a number of percolations connects top and bot- 
tom. Further straining the network only results in stretching 
the (fixed number of) percolations. In that case one would 
expect the stiffness to converge to a fixed 'steady state' value. 
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